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ABSTRACT 

Unified dark matter models are appealing in that they describe the dark sector in 
terms of a single component. They however face problems when attempting to account 
for structure formation: in the linear regime, density fluctuations can become Jeans 
stable and oscillate rather than collapse, though it is possible that this difficulty may 
be circumvented by invoking nonlinear clustering. Here we examine the behaviour 
in the fully nonlinear regime, of collapsed objects that should mimic standard dark 
matter haloes. It is shown that the pressure gradient associated with the unified dark 
matter fluid should be significant in the outer parts of galaxies and clusters, and its 
effects obervable. In this case, no flat or falling rotation curve is possible for any 
(barotropic) equation of state with associated sound speed decreasing with density (a 
necessary condition if the fluid is to behave as pressureless matter at high density). 
The associated density profile is therefore also incompatible with that inferred in the 
outer part of clusters. For the prototypical case of the generalised Chaplygin gas, it is 
shown that this limits the values of the equation of state index a that are compatible 
with observations to a 0.0001 or a > 2. This is in line from what is deduced from 
linear analysis. More generally, from the expected properties of dark matter haloes, 
constraints on the sound speed are derived. For the particular case of the generalised 
Chaplygin gas, this further constrains the index to o < 10“® or a > 6.7. For a 
unified dark matter fluid to mimic dark halo properties, therefore, it needs to have 
an equation of state such that the pressure gradients are either minimal or which 
decrease fast enough so as to be negligible at densities characteristic of the outer parts 
of haloes. 

Key words: dark energy - dark matter - galaxies: haloes - galaxies: clusters: general 
- cosmology: miscellaneous 


1 INTRODUCTION 


of dark energy, particularly perplexing is the unexpectedly 
small value and its relatively recent domination of the energy 
budget of our universe. These issues have triggered much 
investigation, seeking a myriad of alternatives to the cos¬ 
mological constant as the driver of late cosmic acceleration 
(e.g., Li et. al. 2011; Yoo & Watanabe 2012). 


Barring significant modifications to known physical law, 
compelling evidence suggests that most of our universe is 
composed of missing dark energy and matter. The dark en¬ 
ergy component is compatible with a cosmological constant 
A, which effectively contributes a negative pressure, equal to 
the energy density in magnitude but reversed sign (p = —p). 
Dark matter on the other hand generally needs to be ’cold’ 
(or at least warm), in the sense of having small enough 
streaming length at matter and radiation equality, so as 
to satisfy structure formation constraints, including Lyman- 
alpha bounds (e.g., Viel et. al. 2013). The cold dark matter 
(CDM) scenario has developed over three decades into a re¬ 
markably successful model for structure formation, despite 
some problems on galactic scales (e.g., Frenk & White 2012 
for a review). 


Dark energy and dark matter need not form separate 
components. In particular, if a single unified dark fluid obeys 
an equation of state such the magnitude of its (negative) 
pressure decreases rapidly enough with density, then it may 
behave as CDM at high densities and much like A at low 
densities. This contention is appealing and plausible in prin¬ 
ciple, given the observation that the negative pressure con¬ 
tribution to the content of the universe only becomes im¬ 
portant in recent cosmological epochs, when the total energy 
density has been significantly diluted by expansion; and that 
the bulk of structure formation, including collapse leading 
to gravitationally bound structures, took place before dark 
energy became dominant. CDM-like hierarchical structure 


Despite the successes of the ACDM cosmology, the ori¬ 
gin of the missing components remains obscure. In the case 
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formation may then proceed naturally, as in ACDM prior 
to A domination. In this picture, the bound structures, thus 
formed, remain at high density, impervious to later cosmo¬ 
logical density dilution, when the unified dark matter gas 
acts as dark energy on large scales, playing the role of A. 

Most extensively studied among unified dark matter 
models has been the generalised Chaplygin gas cosmology 
(Kamenshchik, Moschella & Pasquier 2001; Bento, Bertom- 
lami & Sen 2002). Here, the relevant fluid is endowed with 
the equation of state Q 


The positive parameters A and a can be adjusted so that 
the associated solutions of the Friedman equations match 
the expansion history of the universe, as revealed by super¬ 
novae type la data. Park et. al. (2010) find a preferred value 
of a = 0.997. This is remarkable, as this original Chaplygin 
gas cosmological model (with a = 1) is particularly well mo¬ 
tivated from a theoretical point of view (Gorini et al. 2006; 
Jackiw 2000). However, this value of a turns out to be in 
tension with the estimated age of the universe, and in catas¬ 
trophic contradiction with probes of structure formation. In¬ 
deed, it has long been known that the non-zero pressure as¬ 
sociated with unified dark matter may significantly stem the 
growth of structure: perturbations can become Jeans stable; 
rather than collapse they oscillate, imprinting unobserved 
signatures on the matter density fluctuation power spectrum 
(Sandvik et al. 2004). Including the baryonic component in 
the analysis of the growth of fluctuations ameliorates the sit¬ 
uation somewhat (Bega et al. 2003). Nevertheless, values of 
a very close to zero or greater than one (when sound speeds 
can become superluminal) are still favoured (Gorini et. al. 
2008a: Fabris et. al. 2008; Fabris, Velten & Zimdahl 2010). 
A final apparent coup de grace comes from combining super¬ 
novae data with matter power spectrum constraints, as well 
as cosmic microwave background and baryon acoustic oscil¬ 
lations data. This seems to conclusively constrain a to very 
small values near zero, with the implication that the ensuing 
model is virtually indistinguishable from ACDM (e.g. Park 
et. al. 2010; Wang et. al. 2013). 

The above constraints were deduced by considering lin¬ 
ear growth of structure from adiabatic perturbations. Non- 
adiabatic perturbations (e.g., Reis et. al. 2003; Zimdahl & 
Fabris 2005; Billie, Tupper & Viollier 2007) and nonlineari¬ 
ties were therefore considered. It is not clear whether ini¬ 
tially non-adiabatic perturbations would remain impervi¬ 
ous to pressure forces into the nonlinear regime, and thus 
whether invoking these solves the problem of structure for¬ 
mation in unified dark matter models. Nonlinear collapse 
(including pressure forces) was investigated by Billie et. al. 
(2004), starting from a homogeneous background and us¬ 
ing a simple collapse model coupled to the Press-Schecter 
formalism. They found that only a small fraction of initial 
density perturbations lead to collapsed structures akin to 
dark matter haloes. Their model is nevertheless open to re¬ 
finement. 


^ Physical units are used here instead of the customary c = 1, as 
a principal concern will be to compare typical speeds in galaxies 
with sound speeds in a unified dark matter fluid. 


Another crucial issue concerns the fact that collapsed 
unified dark matter structures (that form bound CDM-like 
haloes), ’decouple’ from the pressure budget, they there¬ 
fore do not contribute to the dark energy content (Billie 
et. al. 2003; Avelino et. al. 2004; Avelino, Bolejko & Lewis 
2014). This component, already collapsed into bound sys¬ 
tems, behaves as pressureless matter, and can thus cluster 
accordingly. Clustering also affects the background evolu¬ 
tion (which is driven by a modified effective pressure), and 
it has been suggested that the resulting effective equation of 
state could allow Chaplygin gas cosmologies to be compati¬ 
ble with available observations for all relevant values of the 
index a, provided the clustering is efficient enough. This is a 
powerful argument, though accurately evaluating the effec¬ 
tive equation of state would require detailed understanding 
of Chaplygin gas hydrodynamics as the fluid forms struc¬ 
tures in a cosmological context, or at least much more de¬ 
tailed modelling of the collapse and clustering process than 
has been hitherto available. 

Nevertheless, whatever the intermediate evolution, if 
the final quasi-equilibrium ’haloes’ are to be indistinguish¬ 
able from their largely successful CDM counterparts, it is 
required that they have similar density profiles in the outer 
regions. For this to happen, the pressure gradients associ¬ 
ated with the Chaplygin gas should either be negligible or 
somehow collude as to result in an empirically viable den¬ 
sity profile, at least up to radii where direct discriminants 
(such as rotation curves of galaxies and X-ray emissions from 
clusters) are available. These requirements may not be easily 
satisfied, as the average density of CDM haloes within their 
virial radii is only ~ 200 times more dense than the average 
density of the universe (by definition) and is smaller in the 
outer regions; pressure gradients could therefore be palpa¬ 
ble in the outer parts of galaxies. On the other hand, the 
possibility that significant pressure gradients could collude 
so as to lead to profiles similar to those of CDM haloes is 
unlikely, given that - as opposed to the case of collisionless 
CDM - the existence of definite equation of state ensures a 
unique solution to the hydrostatic equilibrium (and associ¬ 
ated density profile), given the boundary conditions. 

In this context, we consider the effect of pressure gradi¬ 
ents on the outskirts of spherical unified dark matter haloes 
in hydrostatic equilibrium, showing that the effect of such 
gradients should be observable if present. In the case of the 
generalised Chaplygin gas, this restricts plausible values of 
a. More generally, we provide constraints on possible val¬ 
ues of the sound speed (and therefore equation of state). 
Section gives a brief introduction to the basic properties 
of Chaplygin gas, then discusses its self-gravitating hydro¬ 
static equilibrium and associated density distribution; more 
general properties of any fluid with sound speed increasing 
with decreasing density are discussed in Appendix B. The 
transition from CDM-like behaviour is tackled in Section |3l 
the resulting outer rotation curves are then compared with 
their observed counterparts. In Section |4] further, and more 
quantitative, constraints on plausible values of the gener¬ 
alised Chaplygin index a are obtained by comparing the 
density profiles of Chaplygin haloes with observationally in¬ 
ferred cluster halo profiles. In Section [5] general constraints 
on the equation of state are derived from the general prop¬ 
erties expected of dark matter haloes. 
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2 THE CHAPLYGIN GAS AND ITS 
SELF-GRAVITATING EQUILIBRIUM 


2.1 Basic properties 

The sound speed associated with the equation of state © 
is given by 


2 A 2 —a —1 /r\\ 

Cs = — = Ac ap . (2) 

dp 

As we will see below A takes values that are of the order 
of the present critical density to the power a + 1. Thus we 
can immediately see that the sound speed will tend to zero 
if Q —>■ 0, for any p, or at high densities (relative to the 
critical density) if a is large enough. This property will be 
important in what follows as it will determine when pressure 
gradients are important in haloes. 

Using the equation of state, along with those govern¬ 
ing the evolution of a homogeneous expanding universe, one 
finds that the density of the homogeneous component in 
such a universe evolves with the scale factor as 


P = 



a3(i +“)) 



( 3 ) 


where B is another positive constant, and the scale factor is 
chosen so that a = 1 at the present cosmological epoch. 
A Chaplygin gas dominated universe will therefore have 
p ~ 1/a® for sufficiently small a, as in the standard mat¬ 
ter dominated case. With a present-epoch energy density 
p{a = 1) = Pc = (A -I- B) i+“ , eventually p —)■ A i+“ , and 
the universe is entirely dark energy dominated. It is thus 
possible to associate the constant A with the dark energy 
content, and the second term in the brackets of equation da 
with a matter contribution. This fixes the constant A to 


A = p/+“(i-r!)„+“) 


l+a\ 


( 4 ) 


where Qm is the present matter contribution to the critical 
density. Alternatively, one may leave A as a free parameter 
to be determined by best-htting models to observations (e.g., 
Wu & Yu 2007), in which case one obtains 0.6pc^“ ^ A < 
0.85pc^“; that is, approximately the same values as obtained 
from equation Q. In general, A < p“^® for a universe that 
is on the verge of being dominated entirely by dark energy. 
In the following we use equation ([4)) with Qm = 0.3 and 


3 


( 5 ) 


to determine A. 

Finally, we note that this determination assumed a ho¬ 
mogeneous Chaplygin gas with no clustering. The effect of 
clustering on hxing the values of the Chaplygin gas param¬ 
eters has not been studied in detail. The simple clustering 
model of Avelino et. al. (2014) predicts that the effective 
equation of state for the homogeneous component of a highly 
clustered system becomes, asa—>■!, ^ = —1 (in units of 

c = 1), which implies that A = p°U^■ At a = 1, their Fig. 1 
suggests that p_ ~ 0.75pc if one starts with 90% clustering 
and with 0 = 1. This means that A ~ O.Op®, which is of the 
same order as the estimates above. 


2.2 Hydrostatic equilibrium 


In the standard (CDM) scenario of structure formation, 
haloes are formed as primordial density fluctuations turn 
around, collapse, and virialise (e.g. Mo, van den Bosch & 
White 2010). During this process the dynamics is collision¬ 
less and so pressure forces are absent. Pressure forces are 
important for the dynamics of the gaseous baryonic com¬ 
ponent, which obeys a conventional equation of state with 
positive pressure growing with density. This means that if 
pressure forces are non-negligible at the initial stage of a col¬ 
lapsing gaseous sphere, they can only increase as this sphere 
shrinks. The situation is reversed if the sphere is made of 
unified dark matter, where the magnitude of the pressure 
decreases with density. In this case, if the pressure forces are 
negligible at the initial stage of the collapse, they will remain 
so as the system shrinks. This property may facilitate the 
collapse, but could also render pressure forces important in 
the outer regions (with small densities) of the final structure 
that ensues from that collapse. 

Before examining the possible effects of pressure gradi¬ 
ents on density distributions of outer parts of haloes, con¬ 
sider first the hydrostatic equilibrium of a sphere of gener¬ 
alised Chaplygin gas. Such an equilibrium can be achieved 
because, although the pressure p is negative, it is more neg¬ 
ative for smaller densities; it is thus an increasing function 
of density, and so its radial derivative is negative when the 
equilibrium involves a negative radial density gradient. Such 
a pressure gradient can thus act to balance self-gravity. 0 For 
velocites and densities expected in dark matter haloes, the 
equation for hydrostatic equilibrium reduces to the standard 
Newtonian form, which equates the pressure gradient to the 
gravitational force at radius r (Appendix 1^: 



(6) 


where M = M(< r) = 4?! p{r)Udr. Using ([T}, one can 
derive 


,dp q-i-2 GM 
aA— = -p ^ 
dr r-^ 


( 7 ) 


-2 

which has a power law solution p ^ r . For small a this 
density profile is similar to the central cusp of an NFW halo 
(Navarro, Frenk & White 1997). And it comes with a rising 
rotation curve for any positive a. 

This is generic; as shown in Appendix no flat or 
falling rotation curves - that is, ones associated with den¬ 
sity distributions falling as p ~ 1/r^ or steeper- can be as¬ 
sociated with pressure supported systems with sound speed 
decreasing with density. Thus, unless pressure gradients are 
negligible throughout the observable regions of galaxies and 
clusters, or the form of the equation of state is modified 
in the presence of significant density gradients, such models 
are conclusively incompatible with observations. Most of the 
following is concerned with examining the conditions under 
which the flrst possibility can be realised. As pointed out in 
Appendix m the second possibility does not appear viable. 


^ Note that in the case of positive pressure decreasing with den¬ 
sity, balancing gravity would require a positive density gradient; 
i.e., density increasing with radius, which is not observed in stellar 
or galactic sysyems. 
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r[kpc] 


Figure 1. Outer density profile of a Milky Way like halo in the 
presence of Chaplygin gas pressure gradients. The halo properties 
are normalised to current estimates of the dark matter density at 
the solar neighborhood. The profile flattens at large radii due to 
pressure forces. 



Figure 2. Outer rotation curves of a Milky Way like galaxy in 
the presence of generalised Chaplygin gas pressure, as function of 
the index a of equation (1). The rotation curves rise due to the 
effect of pressure gradients. 


d(aV) _ _ GM 
dr ^ 


This is the same as equation with ^ replaced by , 
where where cr = (r(r) is to be associated with the one di¬ 
mensional velocity dispersion supporting a collisionless equi¬ 
librium. However, whereas equation m has a unique solu¬ 
tion, given the equation of state and boundary conditions, 
equation ([§1 supports an infinity of density and velocity dis¬ 
persions pairs that solve it. It is indeed remarkable that out 
of all possible solutions, cosmological halo profiles fall into a 
highly constrained set, reasonably well approximated by the 
NFW fitting formula. This is not expected to be reproduced 
by the unique solution corresponding to the Chaplygin gas 
at all radii (as we saw in the previous section). 

In general, both the effective pressure, born of the veloc¬ 
ity dispersion of fluid elements, and the pressure emanating 
from the equation of state may contribute, so that 


^ d(aV) _ _ GM 
dr dr ^ 


( 9 ) 


As long as ^ ^ equation reduces to the Jeans 

equation with associated CDM-like behaviour. However, if 
^ oil® expects the solution to tend towards that 

of the hydrostatic equilibrium of the previous section, with 
its rising rotation curve. In the regions of flat rotation curve 
haloes are isothermal, so that a is nearly constant. In such 
regions, the condition for collisionless behaviour becomes 


dp 2 .. 2 

— =Cs-<u ; 

dp 


( 10 ) 


meaning that, for a given self gravitating mass density dis¬ 
tribution, the speed of sound of the unified dark matter fluid 
making up the halo has to be much less than the equilib¬ 
rium velocity dispersion that could support the density pro¬ 
file against gravity. This condition should also hold approxi¬ 
mately in the outer parts of NFW haloes where the rotation 
curves are slowly falling with radius. 

Flat rotation curves are are associated with density pro¬ 
files falling as ~ 1/^, and in a collsionless systems these are 
associated with constant velocity dispersion a. In a system 
made of unified dark matter, additional pressure forces may 
contribute, corrupting the ~ 1/r^ ’isothermal’ profile; since, 
as we saw in the presbious section, when this pressure dom¬ 
inates, the density profiles fall less steeply with radius. To 
check if fiat rotation curves remain uncorrupted over observ¬ 
able radii, one can solve eqnation with constant cr and 
look for deviations from an isothermal profile, born of non- 
negligible pressure gradient 0. If cr is constant and p is solely 
a function of density, one can differentiate with respect 
to r, while making use of the relation dM/dr = Aixr^p, to 
obtain 


d^ p X / dp\^ 2 dp 
dj -2 Y \ dr) ^ r dr 



( 11 ) 


3 TRANSITION FROM CDM-LIKE 
CONDITIONS TO HYDROSTATIC 
PRESSURE SUPPORT 

A unified dark fluid that behaves as pressureless dark matter 
on galactic scales would be governed by the Jeans equations. 
That is, for spherical equilibrium structures with isotropic 
velocity dispersion 


® Little is known of the dynamics of the transition from colli¬ 
sionless to hydrodynamic behaviour in a Chaplygin gas, it seems 
likely however that the velocity dispersion would decrease con¬ 
siderably as pressure forces become significant (as in standard 
hydrodynamics transition). Hence the assumption of constant a 
may actually underestimate the density (and overestimate the ra¬ 
dius) at which pressure gradients should have significant effect on 
the rotation curves. 


© 2002 RAS, MNRAS 000, [T]-?? 

















Unified dark matter: constraints from galaxies and clusters 5 


Here X = “ p ^ ^ ~ generalised 

Chaplygin gas X = ^ and Y = aAU + 


This procedure is similar to that used in obtaining the 
Lane-Embden equation for stellar polytropes, though we 
skip here the step of transforming to dimnensionless vari¬ 
ables here, keeping the physical quantities which, for illus¬ 
tration, will be directly compared with the specific case of 
a Milky Way like galaxy counterparts. Equation GD) was 
thus integrated while fixing the parameters to values rele¬ 
vant to a Milky Way like galaxy. In the absence of pres¬ 
sure forces, the halo is assumed to be a singular isother¬ 
mal sphere, with p(r) = ^ chosen so that the 

circular rotation speed (vf = 2<j^) is 220 km/s. The corre¬ 
sponding density at solar radius is consistent with the local 
dark matter density inferred from observations and expected 
from standard CDM halo profiles (e.g., Salucci et. al 2010; 
Bovy & Tremaine 2012); if slightly on the higher end, as the 
asymptotic rotation speed is assumed here to arise entirely 
from the gravity of the dark matter component (figures 16 
and 17 of Xu et. al. 2008, for example, suggest that a baryon 
contribution to the rotation speed of 15 — 20% at 60 kpc, 
which translates in a factor of about 1.4 in density). At 
r ~ 10 kpc, the pressure forces are negligible, so that one 
can use the density for the isothermal sphere (and the asso¬ 
ciated radial derivative) to define the initial value problem 
and integrate equation (HU outwards. Note that since equa¬ 
tion HU is local, the fact that a singular isothermal sphere 
may not correctly reproduce the dark matter mass distribu¬ 
tion at very small or large radii is not important, as long we 
are interested in a region where a collisionless dark matter 
density should fall off as 1/r^. 


The results are shown in Fig. [T]for the Chaplygin gas 
(a = 1). Condition (HOI) starts to be violated at r ~ 50 kpc; 
accordingly, the density profile starts deviating from isother¬ 
mal as expected. The corresponding rotation curves are 
shown in Fig. [5] for several values of a. As can be seen, 
unless a is very small (associated with a ACDM cosmology) 
or greater than one (and thus can correspond to superlumi¬ 
nal sound speed), the rotation curve sharply deviates from 
flatness, rising steeply at radii significantly smaller than the 
expected Milky Way halo virial radius (~ 200 kpc). This 
is a consequence of the form of the sound speed as defined 
in equation (El); the multiplicative factor a ensures a small 
sound speed is relatively small at all densities for small a, 
while the sound speed decreases steeply with density for 
large enough a. 


For the Milky way there is no evidence that the rotation 
curve takes to rising at such radii. Indeed, if anything, all 
the evidence suggests the contrary (Xue et. al. 2008; Bhat- 
tacharjee, Chaudhury & Kundu 2014). There are also many 
galaxies where direct measurements of the rotation curve ex¬ 
clude the possibility, while still exhibiting Cj/cr^ > 1. Some 
examples are shown on Table 1. (As in the case of the Milky 
Way above, the rotation curves are assumed to be entirely 
dominated by dark matter at the largest radii probed. We 
thus use = |iic, where Uc is the observed circular ve¬ 
locity at the maximum radius probed. The sound speed is 
calculated from equation (2] using the expected dark fluid 
density of the isothermal sphere associated with cr. Includ- 


Object 

-^^max 

cr 

Cs 


NGC 3741 

7 

32 

51 

2.5 

Malin 1 

88 

151 

352 

5.4 

U1230 

35 

72 

244 

11.5 

U5005 

28 

70 

166 

5.6 

U5750 

22 

56 

161 

8.3 

NGC5055 

48 

120 

168 

2 

UGC 11852 

90 

117 

626 

28.6 

N3198 

44 

106 

181 

2.9 

D154 

10 

29 

125 

18.6 


Table 1. Examples of the velocity dispersions cr, required to sup¬ 
port an isothermal halo profile with fiat rotation curve, exceeded 
by the Cahplygin gas (a = 1) sound speed Cs at maximum ra¬ 
dius, itmax, probed observationally. Data for NGC 3741 is taken 
from Gentile et. al. (2006); Malin 1, from belli, Fraternali & San- 
cisi (2010); U1230 and U5005 from de Blok & McGaugh (1997); 
U5750 from McGaugh, Rubin & de Blok (2001); NGC5055, from 
de Blok et. al. (2008); UGC 11852, from Noordermeer et al.(2007); 
N3198 and D154, from Bottema & Pestana (2015) 


ing a baryonic contribution would decrease the needed dark 
sector density, decreasing cr and further increasing Cs). 


4 CLUSTERS OF CHAPLYGIN GAS? 

The radially extended hot gas in clusters of galaxies allows 
for the probing of the gravitational potential up to large 
radii. When the density profile is decomposed, the dark halo 
density contribution can be reasonably fit using the standard 
NFW profile. Given this, the pressure forces associated with 
Chaplygin gas halo need to be negligible, contributing in¬ 
significantly to the support against self gravity out to the 
maximum radii probed (so that the configuration behaves 
as a collisionless equilibrium that can be fit by NFW). This 
section is concerned with checking whether this is true for a 
sample of galaxy clusters. This will be done by using NFW 
fits inferred from observations, then calculating the pressure 
force associated with the NFW density at the maximal ob¬ 
served radius, then comparing it with the force of self gravity 
of the inferred halo at that radius. Note that the NFW is 
employed here as a useful parametrisation (as it has been 
extensively used in fitting cluster profiles); the crucial point 
being that it fits steeply decreasing outer density profiles of 
observed clusters, while the profiles associated with a uni¬ 
fied dark fluid with significant pressure gradient do not (as 
they fall off too slowly with radius as seen in the previous 
sections). Any profile that fits the observed mass distribu¬ 
tion equally well would give similar constraints, within the 
errors associated with the particular empirical fit, as those 
deduced below. 

The NFW profile is given by 

Pnfw = (12) 

Tg ' ~ Vs ' 

where po is a characteristic density and rs the radial scale 
length. For a halo concentration parameter Cvir = 

47rr|F(cvu)’ ^ ^ 

where the function in the denominator is given by 
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Figure 3. Ratio of the expected generalised chaplygin gas pres¬ 
sure forces to gravitational forces (equation ^Jl for NFW density 
fits recovered by Ettorri et. al (2010) for sample of galaxy clus¬ 
ters. The ratios are plotted at the maximum radii probed and for 
various values of the index ol of equation 1. 


F{x) = ln(l +a;) - . 

1 + a; 


(14) 


Defining x = r/xs, the pressure gradient of a generalised 
Chaplygin gas with this density distribution is given by 

Fp = ^ = —4lc^—Po “(1 + ‘ix)x°‘~^{l + (15) 

or Ts 

while the self gravity term for such a configuration takes the 
form 


Fg 



-‘iTxGplra 


F{x) 

X^{1 + x)^ 


(16) 


In this context, we define the ratio rj = ^ At large 
radiii, Pnfw ~ 1/r^ - 1/r®, which cannot be maintained 
if the pressure gradient is non-negligible, as we saw in the 
previous sections. Therefore, it is necessary that, at the max¬ 
imum radius probed by the observations, 

aAE X^+°‘il+XmA+^°‘il + 3Xm) ,, , 

^ 4^Gp2+“ EsF{Xm) ^ ^ ’ 

where Xm = rmax/rs, rmax being the maximal radius probed 
by observations. We now wish to test where this is the case 
for actual clusters of galaxies. For that purpose we use the 
results of Ettori et. al. (2010), averaging over the two meth¬ 
ods they employed for deducing the NFW parameter and 
setting Xm = Rsp/xs and Xm = R^ap/xs, as given in their 
Table 2. 

The results are shown in Fig. [S] As can be seen, unless 
the index a is very small (< 0.0001) or large enough (a > 2), 
the ratio rj can be of order one or larger. This is agreement 
with the results of the previous section, and also of linear 
analysis. As opposed to the linear case, however, the degree 
of clustering of the large scale Chaplygin gas does not appear 
to be important, since the inferred value of A does not seem 
to vary much with clustering (c.f. Section 2). 

As will be seen in Fig. [S] for each value of a there is 
significant scatter in xj. This includes intrinsic scatter in the 



Figure 4. Ratio of the expected chaplygin gas (a = l)pressure 
forces to gravitational forces (equation [nil for NFW density fits 
recovered by Ettorri et. al (2010), as function of the maximum 
radii probed normalised to the scale radius Xs. The error bars 
represent the maximum possible errors resulting from uncertaini- 
ties in both the measured values of both Xs and concentration 
parameter Cvir- 


values of the inferred halo parameters; significant errors en¬ 
tering in the estimation of those parameter; as well as sys- 
tematics due to the variation of the maximum radius probed 
by the observations. We examine in more detail the impor¬ 
tant case of a = 1 in Fig. |4] and [5l From those figures one 
infers that there is systematic increase in p for clusters that 
are observed out to larger radii. It is therefore safe to as¬ 
sume that- despite the large errors in p, stemming from the 
collusion of errors in estimates of NFW halo scale lengths 
and concentration parameters- the constraints on a are in 
fact rendered conservative by the limitations of maximum 
radii probed. Probing to larger radii is likely to systemati¬ 
cally increase the deduced values of p, further constraining 
the model. In the next section we show that this is the case, 
as tighter constraints are inferred from NFW halo properties 
at the virial radius. 
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Figure 5. Same as in Fig. |4]but with maximum radii normalised 
to the virial radius rvir = cVs ■ 


5 UNIFIED DARK MATTER AND HALO 

PARAMETERS 

Direct observations of the internal dynamics of galaxies and 
clusters usually fall well short of the expected virial radii 
of the haloes believed to envelop these systems. Neverthe¬ 
less, weak lensing fields are also well fit by assuming that 
collapsed structures follow NFW-like density profiles (e.g., 
Mandelbaum 2006; Oguri & Hamana 2011); if anything the 
fits are in fact better when these are truncated beyond the 
virial radius, so as to fall faster than the 1/U behaviour 
associated with an outer NFW profile. This is again incom¬ 
patible with the slowly falling density profiles of pressure 
supported unified dark matter fluids. In this section we ob¬ 
tain further constraints, based on properties of NFW haloes 
derived from simulations, by requiring that the unified dark 
fluid pressure gradient be small up to the virial radius. 

In the most general terms, in order for a unified dark 
matter fluid to behave as CDM, one requires that ^ <C 
^ . For an NFW halo and from r > rs out to the virial 
radius, is a slowly varying function of r and so the above 


condition reduces to that of equation (nni, meaning that the 
speed of sound must be much less that the isotropic equi¬ 
librium velocity dispersion that could support the density 
profile against gravity. At large radii, the velocity dispersion 
of an NFW halo tends to an asymptotic value such that 
(T^ = (e.g., El-Zant 2006, Appendix B), so that 


2 2^F(x) 

<j =TTporsG -. (18) 

X 

We can eliminate by using (1 1311 . In addition Cvir and 
Mvir are found to correlate in numerical simulations, namely 
through power law relations of the form 



(19) 


The logarithmic slope b is found to be of the order of 0.1, 
and the parameters Ah and Bh take values such that c 10 
for a halo of 10*^^ solar masses. For example Klypin, Trujillo- 
Gomez, & Primack (2011) obtain Ah = 9.6 and h — 0.075 
for Bh ~ . (There is also an apparent upturn in the 

mass concentration relation at high masses, though it has 
been argued that this is due to the inclusion of non-relaxed 
systems; Ludlow et. al. 2012). 

Using these relations one finds that 


= UM ^ 

\po 

where 
_ 2(1 -f 36) 
96 ’ 


F(x) 


xF(cv 




e-2/3 


( 20 ) 

( 21 ) 

( 22 ) 


and A is the halo overdensity within the virial radius, such 
that (at z = 0) Mvir = and cJir = 3po ■ 

For A = 200, and values of Ah, Bh and 6 mentioned above, 
is of order 10^'‘(km/s)^. 

At the virial radius 


2 


(7 


CvirF(Cvir)'“l ' 



(23) 


This is a general constraint that the equation of state of any 
unified dark matter fluid should obey. That is, the sound 
speed should either be small everywhere or decrease suffi¬ 
ciently with density so that the above relation is satisfied at 
the virial radius. For example,for the generalised Chaplygin 
gas with an NFW profile the sound speed is 


2 

Cs 


2 


ac 



Q-l-l 

3(0+1) 
Cvir ’ 


(24) 


(where we have assumed large enough Cvir so that Cvir -b 1 ^ 
Cvir and A ~ pf). For a generalised Chaplygin gas to behave 
as pressureless dark matter with an NFW profile up to the 
virial radius thus requires 


(72 


3F(Cvir) 


Q; + 2-e 


2^1 
'-'VI r -*-5 


(25) 


where we have used po/Pc = ^ > . For 6 = 0.1, e = 2.9, 

and for a Milky Way like galaxy Cvir ~ 10, so that a needs 
to be either very small or significantly larger than — (2 — e) 
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for the above inequality to be satisfied. For example, for Cvir 
one needs either a < 2 x 10~® or a > 3.8 for the ratio to be 
less than one;, for Cvir = 30 the index is further constrained 
to a < 10“® or a > 6.7. 


6 CONCLUSION 

In unified dark matter models both dark energy and dark 
matter are reproduced by one component, which can behave 
as pressureless dark matter at high density and as negative 
pressure dark energy at low densities. In a cosmological con¬ 
text, past the radiation era, this implies a phase of matter 
domination followed by dark energy driven expansion. 

In terms of halo structure, the unified dark matter sce¬ 
nario also has its appeal, in principle: bound object, if they 
form, should have density high enough so that the dark fluid 
they are made of behaves as pressureless dark matter, of the 
sort that the haloes harbouring galaxies are believed to be 
made of. Nevertheless, the densities of dark matter haloes 
is not exceedingly large; the average mass within the virial 
radius is by definition only ~ 200 times larger than the am¬ 
bient density of the universe, it is smaller still in the outer 
regions. One may therefore ask if the pressure forces associ¬ 
ated with the unified dark matter can have significant effect 
on the structure of such haloes at observable radii, and what 
that effect, if present, may be. 

In the prototypical case of the generalised Chaplygin 
gas it was shown that - in order to reproduce the observed 
outer rotation curves in galaxies and inferred outer density 
profiles in clusters - the index a of equation 0 should be 
constrained so that either q > 2 or a < 0.0001. Cosmological 
models with small a are quite close to ACDM. The sound 
speed in large a models can be superluminal, though it has 
been argued that this does not necessarily violate causality 
(Gorini et. al. 2008a). Large a models are also in principle 
observationally distinguishable from ACDM, as they have 
the interesting imprint of describing a background evolution 
that transitions rather abruptly from a matter dominated 
to a de Stter phase (Piatella 2010). 

The bounds derived here are similar to those obtained 
from studying the linear clustering regime (e.g., Gorini et. 
al. 2008), though in that latter case it appears that the con¬ 
straints may be circumvented by invoking the effect of non¬ 
linear clustering on the effective equation of state (Avelino 
et. al. 2014). This effect has not been studied in detail, 
though the simplest model (described in the aforementioned 
reference) does not imply modifications that would affect the 
results presented here significantly (cf. Section 2.1). 

The discrepancies, leading to the aforementioned con¬ 
straints, are generic; they stem from the observational re¬ 
quirement that the outer density profile of galaxies and clus¬ 
ters fall off at least as ~ 1 jE , while those unified dark mat¬ 
ter systems have significantly shallower profiles when the 
pressure dominates. Indeed, for a generalised Ghaplygin gas 
configuration with significant pressure gradient, the density 
falls as 1/r or slower (Section 2.2). In general, it can be 
shown that, for any dark fluid with barotropic equation of 
state and sound speed decreasing with density, profiles as¬ 
sociated with flat or falling rotation curves - and therefore 
also a density distribution characteristic of the outer parts 
of clusters - cannot be supported when the unified dark fluid 


pressure forces are dominant (Appendix B). That the sound 
speed decreases with increasing density is a necessary con¬ 
dition if the fluid is to behave as collsionless dark matter at 
high density. 

Requiring that unified dark matter haloes conform to 
NFW profiles up to the virial radius leads to further, tighter, 
constraints on the sound speed and associated equation of 
state. Although the range of empirical dynamical evidence 
prsented here mostly falls well short of the expected virial 
radii of the haloes of the systems examined, weak leasing sig¬ 
nals, which probe the density distribution to mauch larger 
radii, are well modelled by invoking NFW-like halo den¬ 
sity distributions (e.g., Mandelbaum 2006; Oguri & Hamana 
2011) . The NFW profile serves as a useful paramtrisation of 
the density profile of the gravitationally bound component 
of the mass distribution when fitting the leasing data; it pre¬ 
dicts an outer density profile falling as ~ 1/r®. if anything, 
better fits to the leasing data are achieved by truncating the 
halo profiles used so that the density falls off far faster than 
NFW beyond the virial radius. This is in stark contradic¬ 
tion with the slowly falling density profiles we find when a 
pressure contribution associated with a unified dark fluid is 
significant in the outer regions of self gravitating configura¬ 
tions. These shallow profiles should persist well beyond the 
virial radius typical of a GDM I halo, in fact they should 
continue until the unified dark matter fluid composing the 
’halo’ reaches pressure equilibrium with the ambient homo¬ 
geneous background representing the dark energy. Such a 
situation would appear to be in clear contradiction with the 
lensing data. To fit that data one needs a rapidly falling 
desnsity distribution around in the outer parts of self grav¬ 
itating configurations. 
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APPENDIX A: THE NEWTONIAN LIMIT IN A 
GENERALISED CHAPLYGIN GAS HALO 


The equation for hydrostatic equilibrium of a relativistic 
spherical system composed of a perfect fluid with density 
p, pressure p and mass enclosed within radius M = M{< r) 
enclosed within radius r is (e.g. Weinberg 1972) 


dp 

dr 


G_ 


(p + p/U) 



2GM 

r (-2 


y\ (Al) 


In the case of generalised Chaplygin gas with equation of 
state of the form ©, the rich variety of behaviour in the 
relativistic self gravitating regime has been described for ex¬ 
ample by Neven, Tupper & Viollier (2006) and by Gorini et. 
al. (2008b; 2009). In this paper we consider the Newtonian 
limit; which, as we will see, should apply to configurations 
with properties akin to those of dark matter haloes. 

In a dark halo the circular velocity Vc. is much smaller 
than the speed of light c, so that 


2GM 

r (.2 



< 1 . 


Also, using equation o and @, one finds that 


(A2) 


-P 

pc^ 



(A3) 


since the density within a halo is always much larger than 
the critical density of the universe. Finally, 


47rr®p dvrr^pc 
Mc2 ~ 


M 


(A4) 


Since 47rr®pc is of the same order of the mass of homoge¬ 
neous distribution with density pc enclosed inside a sphere 
of radius r, this is much smaller than the actual mass M (r) 
inside such a sphere. The term in brackets is also small (as 
above), so that hand side is also much smaller than one. 

Under these conditions the condition for hydrostatic 
equilibrium reduces to the usual Newtonian equation 


dp__ GM{< r) 
dr ^ r^ 


(A5) 


APPENDIX B: THE ABSENCE OF FLAT OR 
FALLING ROTATION CURVES 

It is possible to obtain flat rotation curves by deriving the 
equations of motion of the Chaplygin gas from a scalar field 
while assuming that the latter possesses space-like gradients 
in the collapsed (self gravitating) phase, which is thus char¬ 
acterised by an accordingly modified equation of state (e.g., 
Armendariz-Picon & Lim 2005). However this comes at a 
hefty price. The gas is no longer an isotropic perfect fluid, 
the corresponding density distribution has poles and is un¬ 
defined in the central regions of the self gravitating system, 
and it general no longer satisfies the positive definite energy 
condition (Diez-Tejedor & Gonzalez-Morales 2013). In this 
paper, as in the studies cited in the previous appendix, it 
is assumed that the equation of state remains unaltered in 
the collapsed phase. It will now be shown that, in this case, 
no flat or falling rotation curves can be obtained for a self 
gravitating system with non-negligible pressure gradient, if 
the sound speed is a decreasing function of density. 
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Cosmology requires that the unified dark fiuid behaves 
as a pressureless matter component prior to dark energy 
domination. If the equation of state remains unmodified this 
implies that, at high enough density, that the unified dark 
matter fluid composing the halo behaves as pressureless dark 
matter at small radii. Nevertheless, pressure forces can be¬ 
come dominant in the outer, low density, regions of the halo. 
As the rest of this paper shows, this will be the case, for ex¬ 
ample, for a generalised Chaplygin gas, unless its index a 
is very small or significantly larger than one. We will now 
show that in such cases it is not possible to sustain flat or 
falling outer rotation curves. 

When the pressure dominates the support against grav¬ 
ity, the circular speed = GM/r) in a halo with pressure 
obeying the equation an equation of state p = p{p) and as¬ 
sociated sound speed is (using equation [6| 


V 


2 

C 



2 dlnp 
d In r 


(Bl) 


Assuming a real sound speed (c^ > 0) therefore implies that 
the density cannot increase with radius. From the above we 
also have 



Now, assume we have a flat rotation curve up to some radius 
ro (so as to fit observations), the associated logarithmic slope 
around ro should then be —2 (so that p ~ 1/r^, M{< r) ^ r 
and is constant). If the rotation curve is to stay flat, or 
fall, the logarithmic slope should be Si —2 for r > ro. In this 
case, the right hand side of equation (O is non-negative. 
On the other hand, for a flat or falling rotation curve the left 
hand side can only be positive if decreases with radius. 
But if the unified dark matter fluid is to behave as CDM for 
high density must decrease with density, and therefore 
increase with radius. Thus it is impossible to have flat or 
falling rotation curves in this case. 
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